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We present a technique which treats, without approximations, the continuum part of the exci- 
tation spectrum in Random Phase Approximation calculations with finite-range interactions. The 
S_l ' interaction used in the Hartree-Fock calculations to generate the single particle basis is also used in 

the Continuum Random Phase Approximation calculations. We show results for the electric dipole 
and quadrupole excitations in the 16 0, 22 0, 24 0, 40 Ca, 48 Ca and 52 Ca nuclei. We compare our 
results with those of the traditional discrete Random Phase Approximation, with the continuum 
independent particle model results and with those obtained by a phenomenological Random Phase 
Approximation approach. We study the relevance of the continuum, of the residual interaction and 
of the self-consistency. We also compare our results with the available total photoabsorption cross 
section data. 

PACS numbers: 21.60. Jz,24.30.Cz,25,20,Dc 
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C . I. INTRODUCTION 

£Nj ■ In the next few years, radioactive ion beams facilities will provide a large amount of data on unstable nuclei. The 
description of the structure of these nuclei is a challenge for the nuclear many-body theories which have been mainly 
tested on stable nuclei. 

The starting point of our description of nuclear systems is the many-body Schrodinger equation with a two-body 
potential built to describe elastic nucleon-nucleon (NN) cross sections and deuteron properties [l, 2]. To this two-body 
potential we add a three-body force whose parameters are chosen to reproduce the 3 H binding energy 3, 4|. Modern 
microscopic calculations, which solve the many-body Schrodinger equation without approximations, describe well the 
structure of nuclear few-body systems as well as that of light nuclei [5, 6]. These results establish the validity of 
the non relativistic description of atomic nuclei. Unfortunately, the computational complexity of these microscopic 
calculations hinders their application to medium and heavy nuclei. 

In recent years, there have been great advances in theories which solve the many-body Schrodinger equation with 
microscopic interactions by doing well controlled approximations 043- The results are very promising but the 
calculations are still computationally rather involved. Probably, the development of these approaches will not be 
rapid enough to cover the requirements for the description of the data which will appear in the near future. We think 
that, for this purpose, effective theories will play a fundamental role. 

Effective theories search for solutions of the many-body Schrodinger equation in a subspace of the full Hilbert space. 
Usually this subspace is chosen to be formed by Slater determinants. This restriction requires the modification of the 
NN interaction in order to reproduce the energy eigenvalues of the microscopic theory. By using Feshbach's projection 
techniques, it is possible to obtain a formal expression which relates microscopic and effective interactions [13j . In the 
common practice, the effective interaction is parametrized, and the values of the parameters are chosen to reproduce 
some experimental data. 

For example, in the Julich approach [lil Il5j|. which we shall call here phenomenological Random Phase Approxi- 
mation (RPA) approach, the single particle (s.p.) wave functions are produced by a Woods-Saxon potential whose 
parameters are chosen to reproduce at best charge radii and s.p. energies near the Fermi level of the nucleus under 
investigation. The values of the interaction parameters are selected to reproduce the excitation energy of some particu- 
larly collective state, for example, in the 208 Pb nucleus, that of the low-lying 3" state. The phenomenological approach 
has been applied in the '80s of the past century to describe, and predict, with success the excitation of the low-lying 
spectrum and of the giant resonances mainly in doubly-magic nuclei. This phenomenological approach, based on the 
Landau-Migdal theory of finite Fermi systems [16], requires the knowledge of a quite a number of observables in order 
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to select the RPA input parameters, i.e. the s.p. basis and the effective interaction. The philosophy of the approach 
requires to change the input for every nucleus considered, therefore, despite of its success, the phenomenological RPA 
approach is not suitable to predict the structure of experimentally unknown nuclei. 

For this purpose, self-consistent RPA approaches are more useful. In these approaches, s.p. wave functions and 
energies are obtained by solving the Hartree-Fock (HF) equations. Since the effective interaction used in HF is also 
used in the RPA calculation this type of calculations is called self-consistent. In the self-consistent approaches the 
values of the parameters of the effective interaction are chosen to reproduce binding energies, and charge radii, of a 
large number of nuclei. This fit produces an effective interaction to be used in all the regions of the nuclear chart, 
including those so far unexplored by the experiment. 

Self-consistent RPA approaches have greater prediction power than phenomenological approaches, but they require 
higher level of accuracy and stability of the calculations. For example, the dimension of the s.p. configuration space, 
beyond a certain size, is not a problem in the phenomenological approach since the effects of the truncation of the s.p. 
basis are taken into account by changing the values of the interaction parameters. This procedure cannot be used any 
more in self-consistent approaches, since the interaction parameters are chosen once forever in HF calculations. This 
drawback of the self-consistent RPA approach is avoided if the full s.p. configuration space is used in the calculation. 
This implies a proper treatment of the continuum part of the s.p. spectrum. 

Self-consistent RPA approaches which consider the continuum have been proposed already in the second half of 
the '70s [l7l - [23j . but they are applicable only if zero-range interactions are used. In this case the continuum RPA 
(CRPA) equations simplify, since direct and exchange matrix elements assume the same expressions, i. e. they are 
proportional. 

Zero-range effective interactions have the great merit of simplifying the calculations. There are however various 
drawbacks in their use, many of them discussed already in Ref. [24] where the Dl parameterization of the finite-range 
Gogny interaction was proposed. Here we would like to mention some other aspects more directly related to the 
present work. In RPA calculations zero-range interactions produce more collectivity than finite-range interactions. 
The difference becomes larger with the increasing value of the momentum transfer [25|. In addition, finite-range 
interactions provide a better description of unnatural parity excitations [TBI . [26l l27l ] . Finally, finite-range interactions 
are more directly comparable with microscopic NN force. 

The reasons presented above motivated our work. We developed a formalism capable to solve the CRPA equations 
with finite-range interactions. The main problem we had to face was the treatment of the exchange terms of the 
interaction matrix elements. In the literature there are few examples of CRPA calculations done with finite-range 
interactions [HI, [29j] , and, to the best of our knowledge, only a single case of self-consistent CRPA calculation [3(| ■ 
The relativistic self-consistent calculations of Ref. [31[ use finite range interactions, but in the so-called Hartree 
approximation where the exchange terms are neglected. Our approach, which will be described in Sect. [XXI follows the 
lines of that proposed by the Jiilich group at the beginning of the '90s [2{|. There are, however, important differences 
due to the different manner of generating the s.p. bases. In the Jiilich case they are produced by using a Woods-Saxon 
potential while in our approach with a HF calculation. 

Details and basic ingredients of the calculations, such as interactions, expansion basis, test of convergence etc., 
are presented in Sect. IIII1 A discussion of a selected set of results is done in Sect. IIV1 We have calculated charge 
conserving excitations in three oxygen isotopes and in three calcium isotopes. We compare our CRPA results with 
those obtained from discrete RPA calculations, with the results of the phenomenological RPA approach and also with 



results obtained by switching off the residual interaction in the RPA calculations. Following the notation of Ref. [32 1 
we shall label independent particle model (IPM) the results of this last type of calculations. In Sect. |V]we summarize 
the main results of our work and we draw our conclusions. 



II. FORMALISM 



The RPA theory describes the excited state of a many-body system as linear combination of particle-hole (ph), and 
hole-particle (hp) excitations. This implies the existence of a s.p. basis which, in our calculations, is generated by 
solving the HF equations: 

H [Mr)] = -^V 2 fe (r) + U(r)<t> k (r) - [ d 3 r'W(r,r')^(r') = e k (f> k (r), (1) 
2m J 

where we have indicated with J£ and <f> k the s.p. hamiltonian and wave function respectively, with IX the Hartree 
potential 

U W = E E / dV^(r')y Q (r,r')^(r'), (2) 

a=l j<F 
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with W the Fock-Dirac potential 



W(r,r') = ^ ^ 0*(r')^(r,r')^(r), (3) 

«=1 j<F 

and with e k the energy eigenvalue of the fc-th s.p. state. In the above expressions, the sums on the s.p. states are 
restricted to those with energy lower than the Fermi energy, i.e. to the hole states. 
In our calculations we consider a two-body NN interaction of the form 

VafaTj) = v a {\r i -r j \) Ofj, with a = 1,2,..., 8, (4) 

where v a are scalar functions of the distance between the two interacting nucleons, and O a indicates the type of 
operator dependence 

Ofj : 1 , r(t) • t(j) , ■ , <r(i) ■ r(i) ■ r(J) , 

S(i,j) , S(i,j)r(i) ■ r(j) , l l3 ■ , ly • r(i) ■ r(j) . (5) 

In the above expression we have indicated with <x the Pauli matrix operator acting on the spin variable, with r the 
analogous operator for the isospin, and with 

S(h3) = 3 , J _ , 2 er(i) -cr(j) (6) 

the usual tensor operator. The terms a — 7, 8 indicate the spin-orbit contributions of the force. We include these last 
terms only in the HF calculations and we consider them in a zero-range approximation as it is done in Ref. [24j , 
We solve the HF equations in a spherical basis, and we express the s.p. wave functions as 

= # Iim (r) = i£ w (r) £ (^^|jm) Xs . (7) 

In the above equation, we have indicated with n,l,j and m the principal quantum number, the orbital angular 
momentum, the total angular momentum and its z-axis projection, respectively. We use r to indicate the distance 
from the center of the nucleus, and fi to indicate the usual angular coordinates of the polar spherical system. The 
symbol ( | ) indicates the Clebsch-Gordan coefficient, Y/ M is the spherical harmonics and Xs the Pauli spinor for the 
spin. We call t the third component of the isospin and we use the convention t = 1/2 for protons and t = —1/2 for 
neutrons. The radial part of the s.p. wave function satisfies the closure relation 
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6 ik Ri(r)R*(r') + ^ R k (r, e k ) R* k (r', e k ) = S(r - r') , (8) 
where we have introduced the symbol 

/> />oo 

% = E + / d£ fc ( 9 ) 



to indicate a sum and an integration on all the s.p. energies above the Fermi surface. In the following, we shall 
indicate explicitly the dependence on the s.p. energies e k which can assume both discrete and continuum values. We 
shall use the index h to indicate all the quantum numbers identifying a hole s.p. state, energy included, while the 
symbol p will indicate all the quantum numbers of the particle state, but the s.p. energy e p . 

A nuclear excited state \v) = \J, II, u>) is characterized by its total angular momentum J, parity II and excitation 
energy u>. In the CRPA theory, the expression of the operator that applied to the ground state generates the excited 
state \v) can be expressed as 



ph 



^ [Xphi e p) a U e p) a h ~ Yph{ e p) a h a p( € p) » ( 10 ) 



4 



where we have indicated with and a the usual particle creation and annihilation operators and with X and Y the 
RPA amplitudes. The CRPA secular equations whose solution provides the values of X and Y can be written as 

(e p - e h - u>) X£ h (ep) + 

E^ [ v ph, P >h>( e p> e P') X p'h'( e P') + u th, P >h'( e P> e P') Y p'h'( e P')] = °> ( n ) 

p'h' e p' 

(e p -e h + u)Y p v h (e p ) + 

E i Htp-h-(t P ^P')Y p v ,M + u^ yh ,(e p ,e p/ )X^ h ,(e p ,)] = 0. (12) 

p'h' e p' 

In the above equations, the interaction terms have been defined as 

v ph,p'h'( e P' e P') — V ph, P 'h'( e P> e P') ~ V ph,p'h'( e P> e P') ' (13) 

and 

Up h ,p> h '(e P ,ep') = {-l) j »' +jh '- J v p T hth , p ,(e p ,e p ,), (14) 

with 

v P Cp'h'( e P' e P') = J dr i r i J dr 2rlRl(ri,e p )Rl,(r2)V p J ^ p " h ,{r ll r 2 )Rh(ri)R P '(r2,e P '), (15) 

v phyh'( e P' e P') = J dr i r i J dr 2rlR* p {ri 1 e p )R* h ,{r 2 )V p J h ^ hl {r l ,r 2 )R p '{r i,e p ,) R h (r 2 ) . (16) 

In an analogous way, according to Eq. (|14[) . we define the corresponding u functions. We have used the following 
definitions for the quantities related to the interaction 

6 

V iffh> ( r i > r 2 ) = ^2v a (n,r 2 ) A J p '™$l, , (17) 

V^ h ,(n,r 2 ) = E V^^)^, (18) 

KCt = ^~ 1 ) jh+jp ' +Kk { % f h , J K } Up3h'K\\V a (Sl)\\j h j P 'K), (19) 

A pk£v = E k { % t K } (3pjh'K\\V a (n)\\jp'3HK) . (20) 



Here for the angular momentum quantum numbers we used the notation a = y/2a + 1. The terms in curly brackets 
are the Racah 6j symbol and the double bars indicate the reduced matrix element as defined by the Wigner-Eckart 
theorem. In our work we adopt the phase conventions of Ref. [33| . 

In the above equations we have factorized the two-body NN interaction (|4]) in a radial part V a (ri,r 2 ), depending 
only on the moduli of the positions of the two-interacting nucleons, and in an angular and operator dependent part 
V a (Ct). We have done this factorization by using the Fourier transformed expression of the NN interaction. 

Our method of solving the CRPA equations consists in reformulating the secular equations (ITT|) and (|T2"j) with new 
unknown functions which do not have an explicit dependence on the continuous particle energy e p . This is the same 
approach adopted in Refs. [H, [29|, but in our case the s.p. wave functions are generated by a HF calculation. For 
this reason, we present here, with some detail, the various steps bringing to the new CRPA secular equations. The 
new unknowns are the channel functions / and g defined as: 

fU r ) = i X; h (e p )R p (r,e p ), (21) 

J Cp 

and 

9ph(r) = i Y p » h *(e p )R p (r,e p ). (22) 
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The first step of this procedure consists in multiplying Eqs. (fTTj) and (fT2|) by R p (r, e p ), the radial part of the particle 
wave function. Considering Eq. (JXJ) , we obtain for the left hand side of Eq. (fTTj) 



(e p -e h -u)R p (r,ep)X£ h (e p ) = "K [R p {r,e p ) X^ h (e p )] - (e h + u)R p {r,e p )X» h (e p ) 



(23) 



The second step of the procedure is to integrate on e p , the particle energy. For the first term in the right hand side 
of the above equation we obtain, using again Eq. (TTJ), 



% [R p { ri e p )X; h {e p )\ 



"£ y2 i R P (r,e p )X^ ph (e p ) + U(r) ^ R p (r, e p ) X; h (e p ) 

f dVW(r,r') R p (r,e p )X; h (e p ) 
-|^V 2 /;,W + U(r)#„(r) J dVW(r,r')/;,(r') 



= M [# h (r)] 



(24) 



In the above equation we have indicated with V 2 the usual Laplace operator where the differential terms related to 
9 and <fr have been already applied to the spherical harmonics providing the correct eigenvalue. Therefore, in our 
writing, we imply that only the derivations on r should be done. 

The operations described above are applied to all the terms of Eqs. (fTT)) and (JT2J) . As example of our calculations, 
we write the contribution of the first interaction term of Eq. (jlll) : 



Rp(r, £p) ^ v p h,p'h' ( e P> e P') X p >h' ( e p) — 

p'h' e p' 

= ^ i? p (r, e p ) £ ^ y dn r 2 | dr 2 r 2 ^(n, e p ) i$ 



^h^fc' (^l , r 2 ) iih (ri ) iZp- (r 2 , ey ) - V^ c ft , (r i , r 2 ) R p > (r i , e p .) R h (r 2 ) 
^ i? p (r,e p ) ^ f dnrj f dr 2 r 2 R* p ( ri , e p ) R* h , (r 2 ) 

^ d p^' Oi ,r 2 )R h (n ) f p , h , (ra ) - V$^ fc , (r, r 2 ) ft , (n ) (r 2 ) 



r J.cxc 



x p 'h'( e p') 



= W dr 2 r 2 ^,(r 2 ) 

p'h' 



V p J C h ,(r,r 2 )R h (r)f», h ,(r 2 ) - V£? h ,(r, r 2 ) / p \,(r) fl fe (r 2 ) 



J,exc 



F(r), 



(25) 



In the above equation we have used the definition f|21[) and the closure relation ([5]) and we have defined the term 



T(r) = - ]T ^i^(r) / dnr^Kn) f dr 2 r 2 ^ ^,(r 2 ) 

e,<e F p'ft' 



^ d p' h> ( r i . r 2 ) R h (n ) f P 'h' (r 2 ) - V^ p 4, (n , r 2 ) (n ) R h (r 2 ) 



where, to simplify the writing, we have dropped the dependence of / on v = (J, II, w). 

We write a new set of CRPA secular equations where the unknowns are the channel functions / and g, 

X [f ph {r)\ - (e h + w)/ ph (r) = -9&(r) + £ 8 ip R 4 {r) f dn r 2 i?*(n) 3 J ph {n) , 
M [/ph(r)]-(c A - c) 3p/! .(r) = -S p 7 ;i (r) + ]T ^(r) | dn r? R* (n ) (n ) , 



(26) 



(27) 
(28) 
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where we have defined 




v P h d p'h> i r ^2)Rh{r) f p > h> {r 2 ) - V^? h , (r, r 2 ) fp> h'(r) Rh (r 2 ] 



+ 



Ip'h' 



(ra) 



^'(^^)i?,(r)i?^(r 2 ) 



(29) 



and is obtained from the above equation by interchanging the / and g channel functions. The relation between 
the U and V symbols is analogous to that of the u and v symbols of Eq. (|14p . The last terms in the right-hand side 
of Eqs. (|2?p and (|28|) are zero if there are no hole states having the same angular momenta / and j of the particle 
state considered. 

We have changed a set of algebraic equations with unknowns depending on the continuous variable e p into a set of 
integro-differential equations with unknowns depending on the distance from the center of coordinates. The solution 
of this problem requires to impose the proper boundary conditions. If the excitation energy ui is above the nucleon 
emission threshold, in some of the ph excitation pairs compatible with the angular momentum and the parity of the 
final state, the particle has positive energy. We call open channels these ph pairs, and closed channels those pairs 
where the particle is in a discrete state. 

After fixing the angular momentum J and the parity II of the excited state, for each value of the excitation energy 
to, we solve Eqs. (|27[) and (128j) a number of times equal to the number of the open channels. Every time we impose a 
different boundary condition, i.e. that the particle is emitted only in a specific channel, which we call elastic channel 
and we label it as po^o- For an open ph channel, we impose that the outgoing asymptotic behaviour of the channel 
function is 



R P0 (r, 



e p ) Sp, Po Sh,h + \H p (e h +uj,r) 



(30) 



where A is a complex normalization constant and H~ (eh + oj, r) is an ingoing Coulomb function or a Hankel function 
in case of a proton or neutron channel, respectively. The s.p. wave function R p is eigenfunction of the HF hamiltonian 
([T]) for positive energy, and is calculated as described in Appendix IBl 

In the case of a closed channel, the asymptotic behaviour is given by a decreasing exponential function 



J ph 



(r — » oo) — > - exp 
r 



2m\e h 



(31) 



as in the case of the channel functions 9^° , 



#°(r^oo) -> - exp 



2m\e h - lj\ 

h 2 



(32) 



We solve the CRPA secular equations (|2"T|) and (|2"8)> by using a procedure similar to that presented in Ref. [29|. 
The channel functions / and g are expanded on a basis of sturmian functions $£J which obey the required boundary 
conditions (|50>(l3"2"j) . 

The sturmian functions $^ are defined as eigenstates of the differential equation [34-38] 



W_ d^_ 

2m dr 2 



to r dr 



h 2 ip{ip + i) 

2m r 2 



(33) 



where to is the particle mass, l p is the orbital quantum number and U p (r) is a complex square well potential of the 
form 



Pg + ijp 1 . if r < a , 
, if r > a , 



(34) 



with f3£ and 7^ real constants. The requirement of continuity of $^ at r = a implies that only a discrete set of values 
of and 7^ should be considered. In this set of solutions, the index /1 is related to the number of nodes of the 
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Sturm-Bessel function $^ in the region < r < a. When the value of the index /x increases by one unity, an additional 
node appears in the wave function at r < a. The definition of the Sturm-Bessel functions given above implies the 
orthogonality relation 



G^+rtf) / drr 2 ^(r)^ p (r) = 5^. 
Jo 



(35) 



Since, in general, the Sturm-Bessel functions are not orthogonal to the wave functions of the s.p. hole states, we 
find more useful to consider a set of orthogonalized functions which we construct as 



W = *£(r) - J2 S *p R W J dr'r' 2 R t (r')^ p (r'), 



(36) 



£i<£F 



where with Si P we indicate that in the sum U = l p and j% = j p . By using this set of orthogonalized sturmian functions 
we express the channel functions fp^ h ° and (& a as 



f v ° h h °{r) = R Po (r,e Po )S PPo S hho + $£+(r) , 



(37) 
(38) 



where the superscripts + and — indicate that the sturmian functions are calculated for e p = e/,. + uj or e p = e/j — u) 
respectively. To simplify the writing we drop the explicit dependence on the open channel label poho of all the c ph 
expansion coefficients. 

We insert the expressions (|37|) and (|38[) in the secular equations (|27|) and (j28|) . and following the steps presented 
in Appendix [AJ we obtain a system of linear equations whose unknowns are the expansion coefficients c p ^ . The new 
CRPA secular equations are 



EE 

H P 'h' 



^ P + ) 



+ E 5 lp (e l -e h -u J )(^; + r\R t )((R i r\^ p + )) 

-(m + rRh'\v p j h %,\R h ^+) - m + rR h '\v p j h ^ h ,\^ + R h )) 
-(((K + r^ P nu£ p ^R h Rh>) - ((^TKnu^iRh'Rh})^, 

= (($; + )*R ho \V p ^ oho \R h R Po (e Po )) - ((% + yR ho \V p ^jR Po (e po )R h ), 



c P 'h' 



1')* 



(39) 



EE 

fl P 'h' 



+ £ M^-^+^)((*r)i^><(^i $ r>; 
-(m-yR h '\v p j h %'\R^) - m-yRh'\v p j c^-R h ) 

{{{^-y^WiC-hWy) - ((^-y^PX^lRh'Rh)) (c%,y 
({^n*R P0 {e P0 )\U-X^ ha \R h R h0 ) - m-yRp (epo)\U^T oho \R ho Rh) ■ 



c P 'h' 



(40) 



In the above expressions, with the bra and ket integration convention we indicate integrations on radial variables 
only. The number of these integrations is given by the number of the functions inserted between the bra and ket 
symbols. For this reason we have inserted the symbol I indicating the identity function. 
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Summarizing, we have converted the CRPA secular equations and (fl~2)) into a set of algebraical equations whose 
unknowns are the expansion coefficients c£ h . These equations have a solution for each value of the excitation energy 
lo above the nucleon emission threshold. 

The solution of the secular CRPA equations provides the channel functions / and g and this allows us to calculate 
the transition matrix elements induced by an operator Tj. If the operator Tj inducing the transition is of one-body 
type of the form 

A 

Tjm(v) = F A r i) OjM^i) 6(Ti - r) , (41) 

i=l 

where we have separated the dependence on the radial and angular parts of the operator, we obtain for the transition 
matrix element the expression 



(J\\Tj\\0) poho = \(ip\\0j\\jh) J &rr*{f;> h °{r)yFj{r)R h {r) 

ph 



+ {-l) J+ S>-i»{j h \\ej\\j p ) / &rr*Rl(r)Fj{r)gH h0 {r) 



(42) 

where with the double bar we indicate the reduced matrix elements of the angular coordinates, as defined in [33] ]. 

In this paper we present results regarding nuclear excitations induced by photons. We consider here only natural 
parity, electric, excitations, and we use the following expression for the operator Tj 

A 

T JM = J2 Zf r( Y JM (Sk) <5(r, - r) , (43) 
i=l 

where Zf is the effective charge 

(£)i[i + „ W ]-(f)i[i-„ W ], »J°-1- 

Zf = 1 (44) 
- [1 + t 3 (i)] , otherwise . 

The second of these expressions is obtained by using an approximation valid for the medium-heavy nuclei we are 
studying [39(. In the above equation A, Z and N are the mass, proton and neutron numbers respectively, and 
Tz{i) = 1 for protons and —1 for neutrons. For a given excitation energy w, and electric transition EJ, we calculate 
the B- value as the incoherent sum on every open channel poho, 

B(uj, EJ : — > J) =J2\(^ J \\Tj\\0) Poho \ 2 . (45) 

Poho 

We obtain the total photoabsorption cross section from the B-value by using the expression 

3 J + 1 e 2 /w\2J-i 



J) = 8V _ p7TW (-) B(u, EJ : — ^ J) , (46) 



where we have indicated with e the elementary charge. 



III. DETAILS OF THE CALCULATIONS 



The formalism developed in the previous section leads to a set of algebraic equations whose unknowns are the expan- 
sion coefficients dt h . The number of coefficients, and therefore the dimensions of the complex matrix to diagonalize, 
is an input of our approach. 

Since the expansion on a basis of sturmian functions is a technical artifact, the solution of the CRPA secular 
equations must be independent of the number of expansion coefficients. We tested the convergence of our results by 
controlling the values of the total photoabsorption cross section in 16 and 40 Ca nuclei. We reached the stability up 
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to the fifth significant figure with 10 expansion coefficients, independently of the multipolarity and of the energy of 
the excitation. 

In our calculations we have used two different parameterizations of the Gogny interaction, the more traditional D1S 
force [4l| and the new DIM force [42[ obtained from a fit to about 2000 nuclear binding energies and 700 charge radii. 
The D1S and DIM forces describe the empirical saturation point of symmetric nuclear matter and reproduce rather 
well the behaviour of the equations of state calculated with microscopic approaches [HI, Hi[ . The situation for pure 
neutron matter is different, because the behaviour of the D1S equation of state at high densities is unphysical. The 
DIM force produces an equation of state which has a plausible behaviour at higher densities, even though it does not 
reproduce the results of modern microscopic calculations. 
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Figure 1: Comparison between the interactions used in our calculations and the microscopic Argonne V8' interaction in the 
scalar channel, vi, spin channel, i>2, isospin channel, V3, and in the spin-isospin channel, V4. The dashed-dotted lines, labelled 
FHNC, have been obtained by multiplying the Argonne V8' terms with the scalar part of the correlation function obtained in 
Correlated Basis Function calculations [1(J |. 



In Fig. [T] we show the D1S and DIM forces in momentum space, and we compare them with the bare Argonne 
V8' interaction [3|. The scalar, v\ (panel (a)), isospin, t>2 (panel (b)), spin, V3 (panel (c)) and spin-isospin, V4 (panel 
(d)), channels are shown. The dashed-dotted lines have been obtained by multiplying the Argonne V8' terms with 
the scalar part of the correlation function obtained in the finite nuclei Correlated Basis Function calculations of Ref. 
[l(J. To be precise we have used the correlation function obtained for the 48 Ca nucleus. In any case, these correlation 
functions are rather similar for all the nuclei considered (see Fig. 21 of Ref. [Uj). The large differences between 
microscopic and effective interactions shown in Fig. [T] indicate that the effective interactions take into account a large 
number of effects explicitly treated in microscopic calculations. The effects originated by the short-range correlations 
are only a limited part of them. It is interesting to notice, in the spin, isospin and spin-isospin channels, the extremely 
large values of the DIM interaction, with respect to those of the other forces. 

In the next section we shall compare self-consistent CRPA results with those of phenomenological calculations. 
The s.p. wave functions of the phenomenological calculations have been obtained by using Woods-Saxon wells, whose 
parameters are given in Ref. [lCf . The residual interaction is a zero-range, density dependent, Landau-Migdal force 
whose parameters are those chosen in Ref. [45| . 

We have investigated nuclei where the hole s.p. levels are fully occupied. This eliminates deformations and minimizes 
pairing effects. 



IV. RESULTS 



A. Oxygen 



We have studied three oxygen isotopes, the doubly magic 16 nucleus, and the 22 and O isotopes. In our model, 
the heavier isotopes are obtained from the 16 core by filling, respectively, the neutron lds/2, and the 2si/2 s.p. levels. 
The ground state properties of these three isotopes, obtained in HF calculations with the D1S and DIM interactions, 
are presented in Table HI and in Figs. [5] and El 
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24q 


16 


22q 


24q 


E/A 


-8.093 


-7.372 


-7.012 


-7.955 


-7.254 


-6.912 


protons 














lsi/a 


-35.37 


-46.43 


-48.64 


-32.74 


-43.38 


-45.65 


1P3/2 


-18.58 


-29.89 


-32.32 


-17.63 


-28.71 


-31.18 


lPl/2 


-12.49 


-23.99 


-25.97 


-11.91 


-23.60 


-25.60 


neutrons 














lSl/2 


-38.61 


-41.02 


-41.11 


-36.00 


-38.32 


-38.46 


1P3/2 


-21.82 


-22.11 


-22.43 


-20.91 


-20.72 


-21.15 


lPl/2 


-15.63 


-18.53 


-17.22 


-15.10 


-17.35 


-16.29 


14/2 




-6.56 


-7.01 




-6.34 


-6.85 


2Sl/2 






-4.13 






-4.09 



Table I: Nuclear binding energies per nucleon, E/A, and s.p. energies of the three oxygen isotopes we have considered, calculated 
within the HF approach by using the D1S and DIM interactions. All the quantities are expressed in MeV. The values of the 
experimental binding energies are -7.976, -7.365 and -7.016 MeV for ie O, 22 and 24 0, respectively. 
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Figure 2: Panel (a): charge density distributions of the 16 nucleus. The shaded area represents the empirical density 
distribution The other lines show the distributions obtained in HF calculations with the D1S and DIM interactions (full 
and dashed lines, respectively) and that obtained with the Woods-Saxon potential of Ref. [10( (dotted line). Panel (b): elastic 
electron scattering cross sections calculated in Distorted Wave Born Approximation by using the charge distributions shown 
in the upper panel, as a function of the scattering angle 6. The empirical charge distribution plotted in the upper panel 
(dashed-dotted line) has been obtained from a fit to the data taken from Ref. [48l. |49|. The numbers in the panel indicate the 
values of the electron energy. 



We show in Table U the binding energies per nucleon, E/A, and the s.p. energies of the three oxygen isotopes. 
The agreement with the experimental binding energies, whose values are given in the table caption, is within the 
few percent. We observe that the DIM interaction is slightly less attractive than the D1S one. In any case, these 
are relatively small differences and we may state that, despite the fact that the two interactions are rather different, 
as we have shown in Fig. [TJ they give very similar descriptions of the binding and s.p. energies of the three nuclei 
considered. The results of the table show that the proton s.p. states become more bound with the increasing of the 
number of neutrons. 

In the panel (a) of Fig. [5]we show the 16 charge distributions obtained in HF calculations where the D1S and DIM 
forces have been used, and we compare them with the empirical charge distribution taken from Ref. (4o| . Our charge 
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Figure 3: Proton p p , neutron p n and matter p m distributions for the three oxygen isotopes we have considered, obtained in HF 
calculations with the D1S and DIM interactions. 



distributions have been obtained by folding the proton distributions with a dipole proton electromagnetic form factor. 
We have verified that more modern, and accurate, form factors do not produce sizable differences in our results. Our 
calculations describe well the empirical charge density on the surface but they underestimate it in the center of the 
nucleus. We show in panel (b) of Fig. [2] the elastic electron scattering cross section calculated in Distorted Wave Born 
Approximation by using the charge distributions shown in the panel (a). We compared our results with the data 
of Ref. [HI |49| . The differences between theoretical and empirical densities show up in the cross sections at large 
momentum transfer. The charge distribution obtained with the Woods-Saxon potential by using the parameters of 
Ref. (To| . does not make a better job in describing the empirical density, which is now overestimated in the nuclear 
center. 

In Fig. [3l we show the proton, neutron, and matter distributions, p p , p n and p m respectively, for the three oxygen 
isotopes we have considered. The matter distribution is the sum of the proton and neutron ones. We do not remark 
relevant differences between the results obtained with the two different interactions. The shapes of the neutron 
distributions show the effects of the filling of the s.p. levels, which are empty in the lighter isotopes. In the 22 
nucleus the level, empty in 16 , is completely occupied. This level gives a contribution mainly on the surface. 

The rms radius of the neutron distribution changes from 2.64 fm in 16 to 3.00 fm in 22 0, for the D1S interaction, 
and from 2.61 fm to 2.97 fm for the DIM interaction. The situation is different in 2 O, where the new s.p. level 
to be occupied is the 2si/ 2 - In this case, the main effect is in the center of the nucleus. The values of the neutron 
distributions rms radii are 3.17 fm and 3.12 fm for the D1S and DIM interaction respectively. It is a relatively small 
change on the neutron distribution surface. 

The proton distributions are interesting since there is no change in the occupation of the s.p. levels in the different 
isotopes, therefore all the differences are produced by the interaction between protons and neutrons. In the panels 
(a) and (d) of Fig. [3] we show the proton distributions of the three oxygen isotopes. We notice that the increase of 
the neutron number produces a change in the interior of the nucleus. The proton lsi/2 s.p. wave function becomes 
wider the heavier is the isotope, and since the normalization is conserved, the value of wave function at the center of 
the nucleus becomes smaller. The rms radius of these distributions changes from 2.20 fm in 16 to 2.45 fm in 24 
almost independently of the interaction used. The relevant lowering of the proton distribution in the nuclear center 
is produced to compensate the increase of the neutron density (see panels (b) and (e)). In the panels (c) and (f) we 
also show the matter distributions and we observe that the differences between the various isotopes are smaller than 
those shown separately by the proton and neutron distributions. 
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So far, we have presented HF results which are related to the ground states properties of the three oxygen isotopes. 
We discuss now the excitation spectra obtained by our CRPA calculations. A first point we have investigated is related 
to the relevance of the proper treatment of the continuum in the self-consistent CRPA calculations. Our study has 
been conducted by comparing our CRPA results with the results of discrete RPA calculations, such as those of Ref. 
(27j . This is the same strategy adopted in Ref. [3(| and, within a relativistic framework, in Ref. 50]. The discrete 
set of s.p. states is obtained by solving the HF equations in a box with bound state boundary conditions. For all 
oxygen isotopes we use a box radius of 12 fm. Larger values of this radius do not change binding and s.p. energies up 
to the fifth significant figure. There is not such a stability for the unbound, e p > 0, s.p. wave functions and energies. 
We have controlled the stability of the RPA results by selecting the maximum value of the particle-hole excitation 
energy, e^ ax , used in the RPA calculation. For a given total angular momentum and parity of the excitation, this 
value determines the number of s.p. states forming the configuration space of the discrete calculation. 

The quantity we have considered for these convergence tests is the centroid energy, which we calculate as 

dujujB(uj,EJ : — > J) 

= -78== • (47) 

/ do; B(lu, EJ : — > J) 

We have studied the convergence for the 1~ and 2 + excitations in all the nuclei we have investigated, and we have 
found that the change from e^ ax =200 MeV to ejj) l ax =250 MeV modifies the value of the centroid energies for less than 
one part on a thousand. All the discrete RPA results we present here have been obtained by using e^ ax =250 MeV. 
Our choice ensures the convergence of discretized RPA calculations done with HF basis which has been generated by 
using a specific value of the box size. 
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Figure 4: Total photo-absorption cross sections calculated with discrete and continuum RPA approaches. The vertical bars 
show the discrete results, and the solid lines those of the CRPA calculations. The excitation multipole is the 1~. The upper 
panels show the results obtained with the D1S interaction, the lower panels with the DIM interaction. 



In Fig. 2] we compare the total photoabsorption cross sections calculated for the excitation of the 1~ resonance in 
the three oxygen isotopes we have studied. The vertical bars show the discrete RPA results and the solid lines those 
of the CRPA calculations. In the upper panels we present the results obtained with the D1S interaction, and in the 
lower panels those obtained with the DIM interaction. 

In the 16 nucleus, the agreement between the results of the two different calculations is rather good. Discrete 
results have their maxima in the same position of those of the continuous solutions. There are peaks around 30, 35 
and 40 MeV which do not have corresponding partners in the CRPA cross sections. The DIM cross sections are 
slightly smaller, indicating, again, that this interaction is less attractive than the D1S force. The situation is more 
complicated in 22 0. Discrete and continuum results have similar structures, but the positions of the peaks are slightly 
different. In any case, the cross sections show a richer structure than in the 16 case. This situation is worsening in 
24 where the peaks of the continuous cross sections do not correspond to those of the discrete calculation. 
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Figure 5: The same as Fig. [4] for the 2 + multipole excitation. 



We show in Fig. [5] analogous results for the excitation of the 2 + resonance. In this case the results of the discrete 
RPA are rather different from those of the CRPA, even in the 16 nucleus. The discrete calculations show clusters of 
peaks not present in the continuous calculations. 



16 Q 22 Q 24 Q 

1" 

RPA-D1S 26.34 24.87 22.82 

CRPA-D1S 27.17 25.11 23.07 

RPA-D1M 26.36 24.72 22.81 

CRPA-D1M 27.23 25.04 23.20 

2+ 

RPA-D1S 28.43 32.91 31.58 

CRPA-D1S 30.47 33.47 32.20 

RPA-D1M 28.19 32.26 30.93 

CRPA-D1M 30.51 33.06 31.06 



Table II: Centroid energies in MeV, Eq. (|47[) . f° r the 1 and 2 + electromagnetic excitations in the three oxygen isotopes we 
have studied, obtained with discrete (RPA) and continuum (CRPA) calculations. 

We show in Table HT1 the centroid energies of the electromagnetic responses obtained in discrete and continuum RPA 
calculations. We have considered for cj m i n the values corresponding to the continuum thresholds. The 1~ resonances 
have been integrated up to w max = 60 MeV, while the 2 + resonances up to w max = 120 MeV. The relative differences 
between these centroid energies are smaller than 2% in the 1~ case. In the case of the 2 + excitation we reach the 
maximum value of 4% relative difference between the DIM results in 16 0. 

While the 1~ response to photon excitations is essentially of pure isovector character, the 2 + response is a combi- 
nation of isoscalar (IS) and isovector (IV) modes. For the 2+ excitation we have separated the IS and IV responses 
by including in the expression of the B(EJ) value of Eq. ([4"5l) the operators 

A 

Tj S m = $>/r JM (^) (48) 

A 

T7m = ^r/W^KsW (49) 
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Figure 6: Isoscalar (IS) and isovector (IV) strengths of the 2 + excitation mode for the oxygen isotopes we have studied. 



with J = 2. We show in Fig. [5] the IS and IV responses for the 2 + excitation of all the three oxygen isotopes we 
have investigated obtained by using the DIS and the DIM interactions. The IS responses are concentrated at lower 
energies and show a sharp peak, while the IV responses are broader and they have less pronounced maxima at higher 
energies. In the 16 nucleus the IS quadrupole resonance has been identified in a scattering processes at a peak 
energy of about 21 MeV to be compared with the peak energies of 21.7 and 21.6 MeV obtained with the DIS 
and DIM forces respectively. The centroid energies of our calculations for the 16 nucleus, calculated for CL> max =100 
MeV are 22.94 and 23.13 MeV for DIS and DIM interactions respectively. If we consider aj max =40 MeV we obtain 
22.12 and 22.24 MeV for these centroid energies. 




Figure 7: Total photoabsorption cross sections for the 1~ (panel (a)) and 2 + (panel (c)) excitations. Experimental data (solid 
squares) are from Ref. In panels (b) and (d), we show the sum rule exhaustion functions, Eq. (|50[) . for the two multipole 

excitations. The grey area indicates the same function calculated for the experimental data. Full and dashed lines show the 
CRPA results obtained with the DIS and DIM interactions, respectively, and the dotted and dashed-dotted lines the IPM 
results obtained with the DIS and DIM HF single particle wave functions, respectively. The horizontal lines indicate the sum 
rule limiting values obtained by using a nuclear matter estimate of the enhancement factor. The values of the TRK sum rule 
is 240 mb MeV. 



Another point we have investigated is related to the effects of the residual interaction in CRPA calculations. In Fig. 
[7] the CRPA results (solid and dashed curves) are compared with the IPM results (dotted and dashed-dotted curves), 
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obtained by switching off the residual interaction in the CRPA calculation, and with the data of Ref. [52[. Since 
more than the 90% of the contribution to the total photoabsorption cross section is given by the 1~ excitation, we 
compare the data with results of this excitation mode. The contribution of the 2 + excitation to the photoabsorption 
cross section is shown in the panel (c) of the figure, and it is one order of magnitude smaller than the contribution of 
the 1 _ mode. In the lower panels we present the sum rule exhaustion functions, 



SR(lu) 



duj'a(uj') 



(50) 



calculated for the cross sections shown in the upper panels. 

The results obtained with the D1S interaction do not show significant differences with respect to those obtained 
with the DIM interaction. Evidently, only the CRPA calculations predict the presence, and also the positions, of the 
resonances. The positions of the peaks are well reproduced for both multipole excitations by the CRPA calculations, 
while the IPM results do not give a good description of the data. The sum rule functions obtained with the IPM 
calculations are smaller than those of the CRPA. The differences are larger for the 1~ excitation than for the 2 + 
excitation. This may be due to the different isospin character of the two excitations, which, as we have already 
pointed out, is mainly IV in the 1~ mode, and a combination of IS and IV in the 2 + mode. In the 1 _ case, the 
energies of the peaks of the CRPA cross sections are larger than those of the IPM results. This because the residual 
interactions are repulsive in the isovector channel. The situation is inverted for the lower energy peaks of the 2 + 
excitation mode. The CRPA cross sections present peaks at lower energies than those of the IPM calculations. This 
indicates that the interactions are attractive in the isoscalar channel. The other, wider, 2 + resonances peaked at 
energies of about 42 MeV, have, instead, IV character, and their energies are slightly greater than those of the IPM 
ones. 

The comparison of the CRPA results with the photoabsorption data emphasizes the well known limitations of the 
RPA description of the giant resonances. The strength is too concentrated in the peak region, and the data show a 
wider energy distribution. The sum rule functions of the panel (b) of Fig. [7] further confirm these deficiencies. Even 
though experimental and CRPA curves seem to have the same limiting values, the CRPA curves saturate much earlier 
than the experimental one. Again, the strength is too concentrated in the resonance region. 

The saturation value given by the Thomas- Reiche-Khun (TRK) sum rule is of 240 mb MeV. The isospin dependence 
of our interactions is responsible for the fact that our calculations saturate at higher values. The calculation of the 
enhancement factor k of the TRK sum rule is rather involved for finite-range interactions I 
of k for the D1S and DIM interactions obtained by a nuclear matter estimate [3(J. We obtain k 
interaction and k—0.50 for the DIM force. These values correspond in 16 to sum rule limiting values of 396 and 
360 mb MeV for the D1S and DIM interaction respectively. As we show in the panel (b) of Fig. [7J these values are 
compatible with the results we obtain with our CRPA calculations. 
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Figure 8: Comparison between the self-consistent CRPA results obtained with the D1S interaction (full lines), the DIM 
interaction (dashed lines) and those obtained with the phenomenological approach of Refs. [S3, HBJ where a Landau-Migdal 
interaction has been used (dotted lines). 



Another issue we have investigated is the capacity of a phenomenological CRPA approach to predict the excitation 
spectra of experimentally unknown nuclei. For this reason we have calculated the 1 _ and 2 + excitations of the three 
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1P3/2 
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-34.91 
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-24.14 


-31.74 
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-37.89 


-37.88 
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-35.60 


-35.71 


lPl/2 


-33.01 


-35.32 


-34.46 
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-22.96 


-22.92 


-23.09 


-22.15 


-21.84 


-22.13 
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-17.54 


-18.27 


-17.40 


-17.03 


-17.69 


-17.89 


14/2 


-15.95 


-17.83 


-18.85 


-15.57 


-17.09 
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Table III: Nuclear binding energies per nucleon, E/A, and s.p. energies of the three calcium isotopes we have considered, 
calculated within the HF approach by using the D1S and DIM interactions. All the quantities are expressed in MeV. The 
values of the experimental binding energies are -8. 551, -8.666 and -8.396 MeV for 40 Ca, 48 Ca and 52 Ca, respectively. 



oxygen isotopes by using the phenomenological approach of Refs. [54J,|55|. In the phenomenological calculations we 
solved the CRPA equations without exchange terms by using a zero-range Landau-Migdal force, whose parameters are 
those of Ref. [45[ . In this approach the s.p. wave functions are generated by solving the one-body Schrodinger equation 
with a Woods-Saxon potential. For all the oxygen isotopes, we used the parameters of the Woods-Saxon potential of 
16 O given in [l(J. In Fig. [3] we compare our self-consistent CRPA results with those of the phenomenological approach 
which arc indicated by the dotted lines. 

The phenomenological results compare rather well with experimental data 52| and with our self-consistent CRPA 
results in the case of the 1~ excitation of the 16 nucleus. The position of the peak coincides with that obtained in 
the self-consistent calculations, and all of them are rather close to the experimental one. We remark, however, that 
the global strength of the phenomenological result is smaller than that produced in the self-consistent approach. The 
reasonable agreement between the CRPA results obtained in the case of the 1~ excitation in 16 is peculiar since all 
the other cases show large differences between phenomenological and self-consistent results. The phenomenological 
calculation predicts the isoscalar 2+ excitation of 16 at 21 MeV energy, but it fails in describing the isovector excita- 
tion at higher energies. The differences between the results of the phenomenological and self-consistent calculations in 
the other two oxygen isotopes are large. The total strengths of the phenomenological cross sections are much smaller 
than those of the self-consistent ones. Great part of the resonance structure of the self-consistent cross sections is 
absent in the phenomenological results. 



B. Calcium 



The same type of investigation done for the oxygen isotopes has been repeated for three calcium isotopes: Ca, 
48 Ca and 52 Ca. The ground state properties of these nuclei are presented in Table Hill and in Figs. [9land[T0l In Table 
IIIII we give the values of the binding and s.p. energies. As in case of oxygen, the agreement with the experimental 
binding energies is within few percents. To be precise, we remark that the D1S results are slightly better than those 
obtained with the DIM force. Also in these calculations the DIM force shows less attraction than the D1S interaction. 
As in the case of oxygen, the proton s.p. levels become more bound with increasing number of neutrons. 

In the upper panels of Fig. [5] we show the charge distributions of 40 Ca and 48 Ca nuclei. We use these charge 
distributions to calculate in Distorted Wave Born Approximation the elastic electron scattering cross sections. These 
cross sections are compared with the experimental data [49|, l56j-|59( in the lower panels of the figure. In the figure we 
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Figure 9: The same as Fig. [2] for the 40 Ca and 48 Ca nuclei. The empirical densities are taken from the compilation of Ref. 



[46l | . The elastic electron scattering data on Ca from Refs. [56H-58|, and those of Ca from Refs. 49, 59], have been rescaled 
to match a unique electron energy of 400 MeV. 



also show, with the dotted lines, the charge densities and the associated cross sections obtained from phenomenological 
calculations done by using the Woods-Saxon potential with the parameters given in Ref. The dashed-dotted 

lines showing the empirical densities are taken from the compilation of Ref. [46} . The results of our HF calculations 
show a better agreement with the data than those of the phenomenological calculations. 
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Figure 10: Proton p p , neutron p n and matter p m distributions for the three calcium isotopes, obtained in HF calculations with 
the D1S and DIM interactions. 



The proton, neutron and matter distributions of the three calcium isotopes are shown in Fig. [TO] The main features 
pointed out in the discussion done for the oxygen isotopes are present also in this case, where there are, however, 
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remarkable differences in the details. In the oxygen isotopes the 2s\/2 state was occupied only by the neutrons in 
the 24 nucleus. In the calcium isotopes we have considered, the 2si/ 2 state is always occupied in both proton and 
neutron cases. The 48 Ca is obtained from the 40 Ca nucleus by filling the neutron I/7/2 state, and the 52 Ca by filling, 
in addition, the neutron 2p 3 / 2 state. In the panel (b) of Fig. [TO] it is shown that the filling of the I/7/2 state increases 
the neutron surface, and leaves practically unmodified the density at the center of the nucleus. The filling of the 
2p 3 / 2 state modifies the neutron density mainly around 2.0 fm, but it produces also a small lowering at the nuclear 
center. These modifications change the proton densities at the nuclear center as it is shown in the panel (a). The 
matter distributions obtained with the DIS interaction, and shown in panel (c), are rather smooth in the nuclear 
interior. The DIM interaction generates in 40 Ca narrow s-waves, and this produces a large proton distribution in the 
nuclear center. The corresponding neutron distribution, panel (e), has a hole in the center, and this compensates the 
peak of the proton distribution and produces a matter distribution rather smooth. The same type of considerations 
can be done also for the distributions of the 52 Ca isotope. In general, we observe that, as in the oxygen case, the 
HF calculations find the optimal matter distributions which is rather smooth, even though the separated proton and 
neutron densities may show some rapid changes. 
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Figure 11: The same as Fig. Ufor the three calcium isotopes. Note that the vertical scale for 40 Ca is half that of 48 Ca and of 
52 Ca. 




The comparison between discrete and continuum RPA results is done in Figs. [TTJand[TS]for the 1 and 2 + multipole 
excitations, respectively. As in the oxygen case, we show the contributions to the total photoabsorption cross sections. 
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The agreement between the results of discrete and continuum RPA calculations is slightly worse than in the oxygen 
case. Both type of calculations produce resonances, however the total strengths predicted by the discrete RPA are 
about one half those of the CRPA. For the calcium isotopes, this result is common to all the nuclei, interactions and 
multipole excitations we have investigated. 



40 Ca 48 Ca 52 Ca 

l - 

RPA-D1S 
CRPA-D1S 
RPA-D1M 
CRPA-D1M 

2+ 



RPA-D1S 


24.66 


27.34 


26.84 


CRPA-D1S 


25.02 


26.89 


26.14 


RPA-D1M 


24.39 


26.70 


26.12 


CRPA-D1M 


27.36 


26.29 


25.49 



22.37 22.27 20.56 

21.89 22.40 20.63 

22.10 21.86 20.24 

22.42 21.77 20.33 



Table IV: Centroid energies in MeV, Eq. (|47[) , for the 1 and 2 + electromagnetic excitations in the various calcium isotopes. 
The discrete RPA results have been obtained by using e")f x =250 MeV. 



We show in Table IIVI the centroid energies of the electromagnetic responses obtained in discrete and continuum 
RPA calculations. As we have done for the oxygen isotopes we have taken for w m ; n the values corresponding to the 
continuum thresholds for the various calcium isotopes. The 1~ resonances have been integrated up to w max = 40 MeV, 
while the 2 + resonances up to w ma x= 100 MeV. The relative differences between these centroid energies are smaller 
than 1% for the 1~ excitations and reach the value of about the 6% for the 2 + states. We have studied separately the 
centroid energies of the IS and IV components of the 2 + excitations. By selecting the same values of Wmin and w max 
we reproduce the energy differences of Ref. [3(| ■ 
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Figure 13: IS and IV strengths of the 2 + excitation mode for the calcium isotopes we have studied. 



The structure of the electromagnetic 2 + strength distributions is analogous to that of the oxygen isotopes. They 
are characterized by a narrow IS peak at lower energies, here around 17-18 MeV, and a much wider IV resonance at 
higher energies. The IS and IV 2 + responses are separately shown in Fig. [T3]for the three calcium isotopes we are 
studying, and for the two interactions we are using. The peak of the IS 2 + resonance has been identified in 40 Ca at 
17.7±0.2 MeV in a scattering processes [60]. This value should be compared with the peak energy of 17.4 MeV in 
CRPA calculations done with both D1S and DIM forces. The presence of an IV 2 + resonance in 40 Ca around 32 



20 



MeV has been indicated in the analysis of radiative neutron capture data [61]. Our calculations produce wide, and 
fragmented, IV resonances. We have peaks around 32 MeV but also around 35 MeV. 




Figure 14: The same as Fig. [7] for 40 Ca . Experimental data (solid squares) are from Ref. [53|. The value of the TRK sum 
rule is 600 mb MeV. 



In Fig. [14] we compare our 40 Ca CRPA results with the total photoabsorption data of Ref. [52| and with the result 
of the IPM calculations. The same observations done for the oxygen case are valid also here. The IPM results are 
unable to describe the experimental cross section. The sum rule functions (I50[) shown in the panels (b) and (d) of the 
figure confirm what we have observed in the oxygen case. The strengths of the IPM results are much smaller than 
those of the CRPA calculations. In the 1~ case, the CRPA sum rule function reaches the empirical value but too 
early with respect to the empirical behaviour. This indicates that the strength is too concentrated in the resonance 
region. 

The value of the traditional TRK sum rule is 600 mb MeV, much lower than the limiting values obtained by our 
CRPA calculations. The nuclear matter estimates of the enhancement factors gives limiting values of the sum rules 
of 990 and 900 mb MeV for the D1S and DIM interactions, respectively. We show in the panel (b) of Fig. [T4l that 
these values are compatible with out CRPA results. 

Also in the case of the calcium isotopes we have compared our self-consistent CRPA results, with those of the 
phcnomenological approach. We obtain results analogous to those shown in the oxygen case. The phenomenological 
results show less strength and structure than the self-consistent ones. 



V. CONCLUSIONS 



In this article we presented a technique to solve the CRPA equations in a self-consistent framework. In our 
calculations we used Gogny-like finite-range interactions containing zero-range density dependent terms. Also the 
spin-orbit term, used only in the HF calculations, is of zero-range type. We have shown results for the A = 16, 22 
and 24 oxygen isotopes, and for A — 40, 48 and 52 calcium isotopes. We have compared our results with the available 
experimental total photoabsorption data. We have studied the need of a proper treatment of the continuum, by 
comparing our results with those obtained by discrete RPA calculations. The need of a self-consistent approach has 
been investigated by making a comparison with the results of phenomenological CRPA calculations. 

We summarize here below the main results of our study. 

- The D1S and DIM forces are very different if compared in the various interaction channels Eq. However, 
they produce very similar results, both in HF and in CRPA calculations. 

- In the HF calculations the global matter distribution, given by the sum of the proton and neutron distributions, 
is modified to obtain the minimization of the total energy of the system. We have shown in Figs. [3] and [10] 
that the HF minimization procedure produces rather smooth matter distributions, even though separately, the 
proton and neutron densities may strongly variate. 
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- With the increasing number of neutrons, the protons s.p. states become more bound, as we have shown in Tabs. 
Uland inTl This effect is not relevant for the charge conserving excitations, those treated in this work, but it may 
have consequences in charge exchange excitations. 

- The IPM calculations of the nuclear responses are unable to provide the proper strength of the multipole 
excitation. This confirms a well known fact that IPM calculations do not predict the presence of giant resonances 
in the nuclear excitation spectrum. 

- The comparison with discrete RPA results shows the need of a correct treatment of the continuum in self- 
consistent calculations. In discrete calculations we have used s.p. configuration spaces large enough to ensure 
the stability of the results in the low-lying states and in the giant resonance excitation regions. Discrete RPA 
calculations can reproduce some bulk properties of the excitation, the centroid energies for example, but they 
fails in the detailed description. Discrete RPA responses show clustering of excited stated that the CRPA 
strength distributions do not have. 

- The phenomenological CRPA results are similar to those of the self-consistent CRPA calculations in 16 O and 
40 Ca . The input parameters of the phenomenological calculations have been chosen to reproduce some ex- 
perimental quantities in these nuclei. On the contrary, when we apply this approach to the other isotopes, we 
found results which are rather different from those obtained with the self-consistent approach. This indicates 
the inadequacy of the phenomenological approach in the study of nuclei lying in experimentally unexplored 
parts of the nuclear isotope chart. 

- Self-consistent CRPA calculations describe rather well the experimental positions of the giant resonance peaks, 
for both the 1~ and the 2 + excitations. On the other hand, the strength distributions are incorrect since they are 
concentrated in the peak region, while the experimental distributions are wider. This is a well known deficiency 
of the RPA description of nuclear giant resonances. There are strong indications that the problem could be 
solved by considering the excitation of two particle-two hole pairs [62h64| . 

The work presented here is the first step of a project aiming to apply a self-consistent computational scheme to 
many other observables and nuclei. The next step of our work will be the study of unnatural parity excitations, and 
for this investigation we shall consider a tensor term in the interaction [(55[ . It has been shown that this term slightly 
affects the ground state properties, [66, 67] but it has more relevant effects on the spectrum of magnetic states |26t l27 
The study of charge-exchange excitations will be a following step. 



Appendix A: Expansion of the CRPA equations on a basis of sturmian functions 

In this appendix we derive Eqs. (l39l) and (l40l) by inserting the expansions (|37| (f38|) of the / and g functions on the 
sturmian functions basis in the CRPA secular equations (|27p and (|28p . We insert the expression ([57)1 in the first term 
of Eq. (|2"T|) and, by using the definition (j3"6")l of the orthogonalized sturmian functions and the fact that R p (r,e p ) is 
an eigenfunction of the s.p. hamiltonian 3-C (see Eq. (fT|)) for the eigenvalue e p = e/, + w, we obtain 



R P0 (r, £p) 5 PP0 S hho + c p t % + ( 
-(e h + u) [R Po (r,e p )S PPo S hho + ^ <+^+(r) 

J2 -*h- u)Ri(r) J dr'r' 2 R t (r')^ p +(r') 



(Al) 



ei<£F 



The sum of the last term is limited to the states below the Fermi surface having the same orbital and total angular 
momentum of the particle state. Using the definition of the Sturm-Bessel functions given in Eq. (1331) . we obtain for 
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the above expression 



E 



ph 



U;(r) - U(r)J $£ + (r) + J dr' r' 2 W(r, /) 

ei<e F J J 

where we have used the fact that, from Eqs. ([T]) and (|33p . we have 

JC[$£+(r)] - pX(r) + e p ] $£+(r) + / dr' r' 2 W(r, r') = - t/£(r) $£+(r) 

Multiplying the above expression by r 2 (r) and integrating on r we obtain: 
drr 2 $;+(r) {:K[/^°(r)] - (e, + w)/$ ft °(r)} = 



E 



y drr 2 $p + (r) 



U(r)$£ + (r) - / dr'r' 2 W(r,r')$£ + (r') 
- E *ip(Q-c/.-w)i2i(r) f dr'r' 2 R t (r')^+(r') 
E c ^ + {^ - <(^ + )*|U|^+> + ((*£+)* I|W|I$£+) 

+ E ^^-e,- W )(^+rii 1 -,)((i? 4 ri^ + >}, 



£i<£F 



(A2) 



(A3) 



(A4) 



where we have used the orthogonality relation ()35j) . The number of the radial integrations is given by the number 
of the functions indicated in the bra and ket symbols. For this reason, in the terms with W we have inserted I to 
indicate the identity function. 

For the right-hand side of Eq. (|27|). using the orthogonality relation ([36| we obtain 



drr' % + (r) 



E fdr>r' 2 R*(r')? J ph (r') 

= J drr 2 %+(r)+ £ S ip R*(r) J dr' r' 2 R t (r') $»(r') 

+ J drr 2 ^+{r) ]T ^i^(r) J dr' r' 2 R*{r')T J ph {r') 
= - | drr^+M^r) . 



(A5) 
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Using now Eqs. ([25]). ((37|) and (|38l) we have 

- J &rr 2 V+{r)T J ph (r) = j dr r 2 % + (r) £ / d/i 



VXp'H'{ry)R h (r)f p , h ,{r') ~ V$£ h ,(rS) # v (r) R h (r') 



U J p C-w{ry)Rh{r)Rh-{r') - U^Q h ,{r,r') R h ,(r) R h (r') 



J.exc / /\ 



= - J &rr 2 <S> v +{r) J dr ' r ' 2 

p' h' 



RUr')V p J Ch'(r,r')R h (r) 



Rpo(r',e Po )5 p , iPo 5 h , M + E c%, d>£+(r') 



Rp (r, e P0 ) V,po h-M + E c %> $ P + ( r ) 



Rh{r') 



U J p C, h ,{r,r')R h {r)R h .{r') - C$£f fc , (r, R h ,{r) R h (r>) 

= -((% + yR h jv p J ^ ho \R hRpo (e PO )) EE^((^ + )*^H^ d ;vi^^ + ) 

p'h> n 

+ ((% + rR ho \V^ ho \R Po (e Po )R h ) + EE s^((*p + )*^^<;7vi^ + ^> 



p'V m 



Putting together Eqs. (|A4[) and (|A7|) . we find a new expression of the CRPA secular equation Eq. (|2T|) 

^ p'h.' i . 

+ E S ip (e i -e h -u;)((% + r\R i )((R i r\^ + )) 

e;<£F 

-(((^+)*i?, i n<, d ; v |i?, 1 ^+) - m + TR h '\vX^\K^R h 

(((% + TKnuX P k>\R^) - {{K + T^pr\u J p p h ,\R hl R h )) (<£v 
= ((^ + )*i?/ io |<^:, |^i? PO ( £po )> - ((^ + )*^ l<C/ l ol^o(^o)^>- 

A similar equation can be obtained from Eq. (|28p for the g channel function. 



c p'h' 



(A6) 



(A7) 



Appendix B: Continuum wave function with HF potential 

We use an expansion on the Sturm-Bessel functions basis to calculate the s.p. wave function for e p > with HF 
mean-field potential. The explicit expression of the differential equation to be solved for the reduced radial part of 
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the wave function u p (r) = r R p {r, e p ) is 

^-j u„(r) + U(r) up(r) - J dr' r'W(r,r')u p (r') = e p u p (r), 



tf_ _ i p (i p + iy 

2m V dr 2 r 



(Bl) 



where U and W have been defined in Eqs. @ and ([3]) respectively. 
We express the solution of the above equation as 



u p (r) 



= Ji P (k p r) +E c p $ pW 



(B2) 



where ji (k p r) is a spherical Bessel function with k p the wave number corresponding to e p . Using the definition (|33p 
of the Sturm-Bessel functions, and their orthogonality relation (|35|) . we obtain the following non homogeneous system 



<W - J drr 2 ^ p {r)U(r)^(r) + j dr r 2 %{r) j dr' r' 2 W(r, r') 

= J drr 2 ^ p {r)U{r)j lp {k p r) - J drr 2 ^ p {r) J dr' r' 2 W(r,r') j lp (k p r') , 
where the unknowns are the expansion coefficients c^. 



(B3) 



Acknowledgments 

We thank I. Sick for providing us with elastic electron scattering data. This work has been partially supported by 
the Spanish Ministerio de Ciencia e Innovation under contracts FPA2009-14091-C02-02 and ACI2009-1007 and by 
the Junta de Andalucia (FQM0220). 



[1] R. B. Wiringa, V. G. J. Stoks, R. Schiavilla, Phys. Rev. C 51 (1995) 38. 
[2] R. Machleidt, Phys. Rev. C 63 (2001) 024001. 

[3] B. S. Pudliner, V. R. Pandharipande, J. Carlson, S. C. Pieper, R. B. Wiringa, Phys. Rev. C 56 (1997) 1720. 

[4] S. C. Pieper, K. Varga, R. B. Wiringa, Phys. Rev. C 66 (2002) 044310. 

[5] S. C. Pieper, R. B. Wiringa, Ann. Rev. Nucl. Part. Sc. 51 (2001) 53. 

[6] H. Kamada, et al., Phys. Rev. C 64 (2001) 044001. 

[7] P. Navratil, J. P. Vary, B. R. Barrett, Phys. Rev. C 62 (2000) 054311. 

[8] D. J. Dean, M. Hjorth- Jensen, Phys. Rev. C 69 (2004) 054320. 

[9] S. Gandoln, F. Pederiva, S. Fantoni, K. E. Schmidt, Phys. Rev. C 73 (2006) 044304. 
[10] F. Arias de Saavedra, C. Bisconti, G. Co', A. Fabrocini, Phys. Rep. 450 (2007) 1. 
[11] S. Quaglioni, P. Navratil, Phys. Rev. C 79 (2009) 044606. 
[12] R. Roth, T. Neff, H. Feldmeier, Prog. Part. Nucl. Phys. 65 (2010) 50. 
[13] P. Ring, P. Schuck, The nuclear many-body problem, Springer, Berlin, 1980. 
[14] J. Speth, E. Werner, W. Wild, Phys. Rep. 33 (1977) 127. 

[15] J. Speth, V. Klemt, J. Wambach, G. E. Brown, Nucl. Phys. A 343 (1980) 382. 

[16] A. Migdal, Theory of finite Fermi systems and applications to atomic nuclei, Interscience, London, 1967. 

[17] G. F. Bertsch, S. F. Tsai, Phys. Rep. 18 (1975) 125. 

[18] S. Shlomo, G. F. Bertsch, Nucl. Phys. A 243 (1975) 507. 

[19] K. F. Liu, N. Van Giai, Phys. Lett. B 55 (1976) 23. 

[20] N. Van Giai, H. Sagawa, Nucl. Phys. A 371 (1981) 1. 

[21] M. Waroquier, et al., Phys. Rep. 148 (1987) 249. 

[22] A. M. Saruis, Phys. Rep. 235 (1993) 57. 

[23] I. Hamamoto, H. Sagawa, X. Z. Zhang, Phys. Rev. C 53 (1996) 765. 
[24] J. Decharge, D. Gogny, Phys. Rev. C 21 (1980) 1568. 

[25] M. Martini, G. Co', M. Anguiano, A. M. Lallena, Phys. Rev. C 75 (2007) 034604. 
[26] G. Co', A. M. Lallena, Nucl. Phys. A 510 (1990) 139. 

[27] V. De Donno, G. Co', C. Maieron, M. Anguiano, A. M. Lallena, M. Moreno- Torres, Phys. Rev. C 79 (2009) 044311. 
[28] T. Udagawa, B. T. Kim, Phys. Rev. C 40 (1989) 2271. 

[29] M. Buballa, S. Drozdz, S. Krewald, J. Speth, Ann. of Phys. 208 (1991) 346. 



25 



[30] H. Nakada, K. Mizuyama, M. Yamagami, M. Matsuo, Nucl. Phys. A 828 (2009) 283. 
[31] J. Piekarewicz, Phys. Rev. C 64 (2001) 024307. 

[32] O. Bohigas, N. Van Giai, D. Vautherin, Phys. Lett. B 102 (1981) 105. 

[33] A. R. Edmonds, Angular momentum in quantum mechanics, Princeton University Press, Princeton, 1957. 
[34] M. Rotenberg, Ann. Phys. 19 (1962) 262. 
[35] S. Weinberg, Phys. Rev. B. 133 (1963) 232. 

[36] S. Weinberg, Lectures on Particles and Field Theory, Brandeis Summer Institute in Theoretical Physics, S. Deser and K. 

W. Ford Eds., Prentice-Hall, Egelwood Cliffs, 1964. 
[37] G. Newton, Scattering theory of Waves and Particles, McGraw-Hill, New York, 1966. 
[38] G. Rawitscher, Phys. Rev. C. 25 (1982) 2196. 

[39] J. M. Eisenberg, W. Greiner, Excitation Mechanisms of the nucleus, North Holland, Amsterdam, 1970. 
[40] J. M. Blatt, V. F. Weisskopf, Theoretical nuclear physics, John Wiley and sons, New York, 1952. 
[41] J. F. Berger, M. Girod, D. Gogny, Comp. Phys. Commun. 63 (1991) 365. 
[42] S. Goriely, S. Hilaire, M. Girod, S. Peru, Phys. Rev. Lett. 102 (2009) 242501. 
[43] A. Akmal, V. R. Pandharipande, D. G. Ravenhall, Phys. Rev. C 58 (1998) 1804. 

[44] S. Gandolfi, A. Y. Illarionov, S. Fantoni, J. C. Miller, F. Pederiva, K. E. Schmidt, Mont. Not. R. Astron. Soc. 404 
(2010) L35. 

[45] G. Co', V. De Donno, C. Maieron, M. Anguiano, A. M. Lallena, Phys. Rev. C 80 (2009) 014308. 
[46] C. W. De Jager, C. De Vries, At. Data Nucl. Data Tables 36 (1987) 495. 
[47] R. Anni, G. Co', P. Pellegrino, Nucl. Phys. A 584 (1995) 35. 
[48] I. Sick, J. McCarthy, Nucl. Phys. A 150 (1970) 631. 
[49] I. Sick, private communication. 

[50] J. Daoutidis, P. Ring, Phys. Rev. C 80 (2009) 024309. 

[51] K. T. Knopfle, G. J. Wagner, H. Breuer, M. Rogge, C. Mayer-Boricke, Phys. Rev. Lett. 35 (1975) 779. 
[52] J. Ahrens, et al., Nucl. Phys. A 251 (1975) 479. 

[53] M. Traini, G. Orlandini, R. Leonardi, Rivista Nuovo Cimento 10 (1987) 1. 
[54] R. de Haro, S. Krewald, J. Speth, Nucl. Phys. A 388 (1982) 265. 
[55] G. Co', S. Krewald, Nucl. Phys. A 433 (1985) 392. 

[56] B. B. P. Sinha, G. Peterson, R. R. Whitney, I. Sick, J. McCarthy, Phys. Rev. C 7 (1973) 1930. 

[57] I. Sick, J. Bellicard, J. Cavedon, B. Frois, M. Huet, P. Leconte, P. Ho, S. Platchkov, Phys. Lett. B 88 (1979) 240. 

[58] J. M. Cavedon, Ph.D. thesis, Universite de Paris-Sud (France), unpublished (1980). 

[59] I. Sick, Model-independent densities of s/d-shell nuclei, unpublished (1975). 

[60] Y.-W. Lui, J. D. Bronson, C. M. Rozsa, D. H. Youngblood, P. Bogucki, U. Garg, Phys. Rev. C 24 (1981) 884. 

[61] I. Bergqvist, R. Zorro, A. Hakansson, A. Lindholm, L. Nilsson, N. Olsson, A. Likar, Nucl. Phys. A 419 (1984) 509. 

[62] S. Drozdz, S. Nishizaki, J. Speth, J. Wambach, Phys. Rep. 197 (1990) 1. 

[63] S. Kamerdzhiev, J. Speth, G. Tertychny, Phys. Rep. 393 (2004) 1. 

[64] D. Gambacurta, M. Grasso, F. Catara, Phys. Rev. C 81 (2010) 054312. 

[65] G. Co', V. De Donno, M. Anguiano, A. M. Lallena. larXiv:1009.336"4l [nucl-th]. 

[66] T. Otsuka, T. Matsuo, D. Abe, Phys. Rev. Lett. 97 (2006) 162501. 

[67] M. Moreno- Torres, M. Grasso, H. Liang, V. de Donno, M. Anguiano, N. Van Giai, Phys. Rev. C 81 (2010) 064327. 



